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A semi-classical model for wobbling motion is presented as an extension to the Bohr-Mottelson 
model of wobbling motion. Using the resultant wobbling potential, a quantum mechanical equation 
is derived for anharmonic wobbling motion. We then attempt to explain the anharmonicity observed 
in the excited bands of two wobbling phonons in the A ~ 160 region. 


Significant progress has been made in a last couple of 
years with regard to the understanding of nuclear wob¬ 
bling motion. This phenomenon was first predicted the¬ 
oretically by Bohr and Mottelson thirty years ago 0, 
but it was not observed until experimental techniques 
advance such as drr gamma-ray detectors. Finally in 
2001, the first experimental report was published on 
the evidence for wobbling motion in ^®^Lu 0]. Sub¬ 
sequently, possible wobbling excitations were reported 
also in 165,167^^ 001. These wobbling bands (of one- 
phonon excitation) were analyzed first by the particle- 
rotor model (PRM) 013 and then by the random phase 
approximation (RPA) [g. 

Evidence for two-phonon bands has also been reported 
in 163.165LU 00. Despite much experimental evidence 
implying that the bands consist of two-phonon exci¬ 
tations, a serious discrepancy with the original Bohr- 
Mottelson model was observed in the experimental en¬ 
ergy spectrum: it shows strong anharmonicity. Cur¬ 
rently, neither PRM nor RPA are successful to explain 
this anharmonicity, but Matsuzaki and Ohtsubo pro¬ 
posed an interesting idea with the tilted-axis cranking 
model (TAG) to account for this anharmonicity 0 . They 
calculated the energy surface with respect to the tilt an¬ 
gles of the total angular momentum vector, and found 
that the curvature of the energy surface around the ori¬ 
gin becomes flatter as angular momentum increases. At 
(and beyond) some critical angular momentum, the lo¬ 
cal minimum at the origin disappears and new minima 
emerge as a manifestation of tilted rotation. They sug¬ 
gested that this “phase transition” from principal-axis 
rotation (PAR) to tilted-axis rotation (TAR) causes the 
anharmonicity, although they did not explicitly demon¬ 
strate a mechanism for the excited bands to acquire the 
anharmonic character. 

The aim of this paper is to present an analytical 
model to explain the anharmonicity observed in the nu¬ 
clear wobbling motion. Two steps are required to derive 
the anharmonic wobbling model proposed in this work. 
First, we consider a semi-classical treatment of the Bohr- 
Mottelson model to derive the wobbling potential energy. 
Then, the re-quantisation is made using the potential, so 
as to obtain the quantum mechanical equation for the 


*e-mail: m.oi@surrey.ac.uk 


anharmonic wobbling motion. 

In the Bohr-Mottelson model, the Hamiltonian reads 


^ = E 



( 1 ) 


This is the Hamiltonian of a triaxial quantum rotor with¬ 
out intrinsic structure. Classically, the dynamics of this 
rotor can be described by the three Euler angles and three 
components of the total angular momentum vector. The 
angular momentum operators A in the Hamiltonian are 
those in the body-fixed frame, so that their commuta¬ 
tion relations are given as where Cijk 

denotes the Levi-Civita symbol. The Hamiltonian com¬ 
mutes with the total angular momentum P = 
but each component A does not commute when the sys¬ 
tem possesses triaxiality {Ji > J 2 > Js, for instance). 

with {i,j,k) being cyclic. In the Bohr-Mottelson model, 
the quantization axis is chosen along the 1-axis. The 
above non-commutation relation implies that the eigen¬ 
value of Ii (denoted in this paper as K) is not a good 
quantum number. Therefore, it is convenient to intro¬ 
duce here the following expression for the wobbling state, 
which is generally described as 

/ 

|wobble;/)= ^ Ck\IK). (3) 

K=-I 


In order to allow the classical analogy, the high-spin 
condition is applied in the Bohr-Mottelson model. Clas¬ 
sically, this condition is expressed as / ~ /i >> 1. Con¬ 
sidering the non-conservation of the K quantum number, 
this classical condition should be interpreted as 

I ^ ih) » 1, (4) 

where the expectation value is taken with respect to the 
state given in Eq. To satisfy this condition, we require 
now 

\Cif » \Ck\^- ( 5 ) 

K^I 

Due to the normalisation condition l^'if P = !> the 
above equation also implies \Ci\^ ~ 1. 
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Now, the Hamiltonian is decomposed into two terms. 

/(/+!) 


H = 


2Ji 




( 6 ) 


where the Dirac constant is set to be unity {Ti = 1). 
The first term is a c-number because / is the quantum 
number (the total angular momentum). The second term 
contains q-numbers, 


i?w = 


1 f 1 


1 


2 \J2 Jl 


1 f 1 


1 


2 \Jz Jl 


P 


(7) 


This approximation is a little drastic from a quantum 
mechanical point of view, not only because the operator 
1 1 is simply replaced by a c-number (if), but also because 
the c-number is a non-conserved if-quantum number in 
this model. However, this approximation can be totally 
justified when we consider the expectation values of the 
relevant quantities. That is, instead of Eqs. and 
we have 


([a, at]) = (/i)//=.l. 

(13) 

(77w) = a;w (n) + ^ 

(14) 


In the Bohr-Mottelson model, bosonic operators of cre¬ 
ation and annihilation are introduced so as to diagonalise 
iiw Two steps are necessary for this aim. First, the cre¬ 
ation and annihilation operators are respectively defined 
as 

o’^ _ ^ 

Second, to eliminate so-called “dangerous terms” (such 
as and aa), the canonical transformation is per¬ 

formed to introduce a new representation of the creation- 
annihilation operators (c, c^). 


- , (15) 

where the expectation value is taken with respect to the 
wobbling state given in Eq. 0 . This “semi-classical” ap¬ 
proximation cannot be compared with the experimental 
energy spectrum at low-lying states directly because the 
number of wobbling phonons can now take any real num¬ 
ber which gives a continuous spectrum. In other words, 
the current model is set back to the classical theory due 
to the semi-classical approximation. 

To re-quantise the system, let us introduce a dynamical 
variable 9 which is defined as 


P = xa^ — ya, 


with 


2 2 1 
X -y =1. 


(9) 


Here x and y are written in terms of the moment of iner¬ 
tia (Ji) and are determined so as to make the dangerous 
terms in the new representation vanish. The commuta¬ 
tion relation is invariant with respect to the canonical 
transformation. The pair of operators (c, c^) therefore 
follows the same commutation relation as the original. 
The resultant Hamiltonian can be written by only c- 
numbers. The corresponding energy spectrum now reads 


Ei 


7(7+1) 

2Ji 


+ Ww 



( 10 ) 


Here n is an eigenvalue of the number operator n = Pc 
and the wobbling excitation energy is given as Ww = 

{I/Jl) ^/{Jl-j2){Jl-Jz)/{j2Jz). 

It looks as if the Hamiltonian is fully quantised, but it 
is the result of the approximation for the bosonic com¬ 
mutation relation. 


9 = cos ^ 



(16) 


This dynamical variable physically means the wobbling 
angle of the total angular momentum vector. If the high- 
spin condition J ~ (Ji) >> 1 is valid, 9 is very small 
{9 « 1), which is the case considered in the original 
wobbling model. The average value for the number op¬ 
erator is calculated as. 


{n + 




(17) 


where the term Re {C 1 C 1 - 2 ) is neglected due to the con¬ 
dition Eq. and K is given as 

I/J2 + I/773-2/J-1 

V'(l/J2-1/J'l)(l/J3-1/Jl)' 

In this way, the mean value of the wobbling Hamiltonian 
is now expressed as 


[a,a/]= ii/I (11) 

which is based on the condition Eq. 0 . The exact form 
of should be, 

Hw = Ww -f . (12) 

Because we can easily prove that n and Ii do not com¬ 
mute with each other, n cannot be treated as a quantum 
number without the approximation in Eq. dJ. 


(77w) (19) 



Because it is expected that the order of k is of 0(1), we 
may justify writing (Hw) as 

(Hw) — ■^Ww(l-cos0) (20) 
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P 

= {IIJ 2 + I/J3 - 2/^1) (1 - COS0) . 

O 

Let us consider the classical motion of 9 by treating this 
energy expectation value as the potential energy for the 
classical wobbling motion. The associated Lagrangian is 
now introduced as 

L{9, 0) = (1 — cos 9). (21) 

The semi-classical oscillator frequency Wsc is introduced 
here as 

= y (1/^2 + I/J3 - 2/Ji). (22) 

This interpretation is possible because under the high- 
spin condition given in Eq. 0 , that is, when the wob¬ 
bling angle is small {9 « 1), the potential is approx¬ 
imated that of a simple harmonic oscillator, that is, 
V{9) ~ . 

Therefore, quantisation for the wobbling motion is car¬ 
ried out through the simple one-dimensional Schrodinger 
equation. 

{-^^+"iWsc(1-cos6I)|t(6») =e;^4'(6»). (23) 

This equation can be considered as an extended wobbling 
model that can handle even the wobbling motion with 
large amplitude (that is, I >> (/i)). It should be noted 
that we did not use the first part of the approximation 
Eq.Q) in the derivation. (The high-spin condition I >> 
1 is used, however.) 

In this model, the mass parameter m is the only 
free parameter. However, using experimental data, 
we can estimate the order of m in the following way. 
Eirst, when 9 << 1, the energy spectrum is given as 
= U}sc{i 2 +^). The one-phonon excitation en¬ 
ergy (AE^'P**™™) observed in ^®^Lu is approximately 
350 keV, which should be equal to Wsc- Using Eq. 
we have 

/M 1 /J 2 + I/J 3 - 2 /^ 1 ) 

m=— - 5 -. (24) 

The anharmonic wobbling potential derived in our 
model, that is, V{9) oc 1 — cos 9, is not strong enough 
to reproduce the anharmonicity observed in experiment. 
According to experiment, the excitation energy from the 
first phonon state to the second is reduced by nearly 
50% in comparison to the excitation energy from the 
yrast to the first phonon state, whereas our anharmonic 
wobbling spectrum gives, at most, 10 % reduction (at 
/ = 30?i), from the result of our numerical calculation 
where the rigid body moment of inertia is used for ^®^Lu 
with (/3, 7 ) = (0.4,—20°). (The degree of triaxiality is 
here measured by 7 |10|. In this study, we employ the 
Hill-Wheeler coordinates, which gives the opposite sign 
convention to the so-called Lund convention.) The mass 


parameter m is positive for the rigid-body moment of in¬ 
ertia with 7 < 35°, which applies to the present case. 
Owing to this fact, we consider only the case that the 
mass parameter is positive (m > 0 ), in the following dis¬ 
cussions. According to the spin-dependence in Eq. ( 1 ^ . 
the mass parameter becomes increased at higher spin. At 
I = 30?i and 40?i, m is calculated as m = 6.7 and 11.9 
(MeV“^), respectively. Consequently, the higher the to¬ 
tal angular momentum, the less anharmonicity. For in¬ 
stance, the reduction is less than 5% at / = 407i. In 
addition, with irrotational-flow moment of inertia , we 
found that the anharmonicity turns out to be much less 
substantial. This is because the moment of inertia is 
proportional to /3^ in the case of irrotational flow, which 
gives rise to one order of magnitude smaller in compari¬ 
son to the rigid body case. As a result, m tends to be one 
order of magnitude larger in accordance with Eq. ( 1 ^ . 

In the rest of this paper, we attempt to explain the 
discrepancy between the above model and observed spec¬ 
trum, which shows very strong anharmonicity. One pos¬ 
sibility for the discrepancy could be the evolution of nu¬ 
clear structure originating from microscopic degrees of 
freedom, which is neglected in the model. For exam¬ 
ple, such an effect can be attributed to the Coriolis force 
in the rotating frame. Not only can the stretching ef¬ 
fect of the nuclear shape cause a change of the moment 
of inertia, but also the quasi-particle excitations due to 
rotational alignment cause a kind of “phase transition” 
from the BCS (or superfluid) phase to the normal fluid 
phase. Therefore, in order to investigate further, we 
might need to go to microscopic theories such as the 
3D-cranked Hartree-Fock-BoTOliubov approach and the 
generator coordinate method |^, which demand massive 
numerical efforts. 

The suggestion by Matsuzaki and Ohtsubo based on 
their microscopic calculations is useful in considering pos¬ 
sible effects originating from intrinsic structure without 
performing demanding numerical calculations. In the fol¬ 
lowing, we focus on a “phase transition” from PAR to 
TAR, as mentioned at the beginning of this paper. In 
order to account for very strong anharmonicity, let us 
consider the following potential, which is similar to the 
energy surface obtained by Matsuzaki and Ohtsubo. 

I4(^) oc 1 — cos 9 — a[l — cos^ 9) . (25) 

FigCl shows the shape of the potential for difference 
choices of the parameter, a. When 9 << I, the above 
potential reduces to Va{9) ~ (1 — 2a)9'^ 12, which implies 
the effective mass m* = (1 — 2a)m. From the simple an¬ 
alytical analysis of the above potential, one can tell that 
the critical point happens at a = 1/2. (See Fignj Below 
the critical point (— 1/2 < a < 1 / 2 ), the classical ground 
state, that is, the minimum of V{9), is found at 6 * = 0. 
This solution corresponds to PAR. Whereas, beyond the 
critical point (a > 1 / 2 ), the classical solutions appear at 
9 = ±cos“^(l/2a), which correspond to TAR. Namely, 
the parameter a has a physical meaning as a control pa¬ 
rameter of a (classical) phase transition between TAR 
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and PAR. (The domain a < —1/2 is not considered in 
this study because an unphysical minimum appears at 
9 = ±7r.) 



FIG. 1: The wobbling potential energy with strong anhar- 
monicity, Va{0), for different control parameter, a. 



FIG. 2: Wobbling spectrum as a function of the control pa¬ 
rameter, a. The rigid-body moment of inertia is employed 
with (/3, 7 ) = (0.4, —20°). 

The physical situation seen beyond the critical point 
(a > 1 / 2 ) is similar to the mechanism of spontaneous 
symmetry breaking in the linear sigma model, or the¬ 
ory in quantum field theory Unlike theory, how¬ 
ever, the cause of the “negative mass” in our case can 
be attributed to more specific physical effects, that is, 
microscopic effects. It is widely known that nuclear mo¬ 
ment of inertia matches neither that of a rigid body, nor 
of an irrotational flow. It takes, in fact, the value be¬ 
tween these two limiting cases The primary reason 
is the presence of the pairing correlation between con¬ 
stituent particles (nucleons). Because the mass parame¬ 
ter m in our rotor model is determined by the moment of 
inertia in such a way presented in Eq. the deviation 
of the realistic value of the moment of inertia from the 


rigid-body or irrotational-flow values can influence the 
magnitude of the mass parameter. It is therefore natural 
that the microscopic effects, such as the pairing correla¬ 
tion, are taken into account through the renormalisation 
of the mass parameter, that is, m*. In other words, the 
microscopic effects are supposed to be included in the 
present model effectively, through the control parame¬ 
ter, a. When the effective mass parameter, m* changes 
from a positive value to negative, one can suppose that 
the spontaneous symmetry breaking occurs as an anal¬ 
ogy with the quantum field theory. A difference from 
the (j)^ theory is that the broken symmetry in our model 
is a discrete symmetry {9 —> —9), so that no Nambu- 
Goldstone mode is created for the restoration of the 
symmetry. Instead, quantum tunnelling plays the role of 
the symmetry restoration. 

To examine how the strong anharmonicity is induced 
through the spontaneous symmetry breaking, numerical 
calculations are carried out for the potential Va{9) by 
varying the control parameter, a. The results are shown 
in Fig.2. Below the critical point (a < 1/2), the en¬ 
ergy difference between the ground and the first excited 
state is of the order of huJsc, which is about 0.35 MeV 
in the present case. The anharmonicity can be seen, but 
they are very small (less than 10 %), as already discussed 
earlier in connection to the original wobbling potential, 
Ea.lEHIl. However, beyond the critical point, the energy 
difference becomes extremely small due to the tunnelling 
effect. As the ground state is bound deeper in the po¬ 
tential, the energy splitting becomes smaller. Finally, 
around a ~ 1 . 1 , the ground and the first excited states 
can be regarded as almost degenerate, which implies an 
onset of the spontaneous symmetry breaking in a quan¬ 
tum mechanical sense. In other words, the quantised 
wobbling motion around 0 = 0° (PAR) goes into a phase 
transition to TAR beyond the critical point (a ~ 1.1). 
Considering that the degeneracy is perfect, the lowest 
three levels (n = 1, 2, and 3) give rise to a very strong 
anharmonicity. In particular, at a = 1.1, the ratio is cal¬ 
culated to be A£’ 2 ^i/AE 3^2 = 335 (keV)/ 167 (keV) 
~ 2, where AEi^j = Ej — Ei. This result agrees well 
with the experimentally observed anharmonicity. 

These ideas are based on the consideration that micro¬ 
scopic degrees of freedom bring further anharmonicity in 
addition to the wobbling potential derived in Eq. (P|l. 
It is thus necessary to check the ideas and models pre¬ 
sented here by using microscopic approaches, which has 
to be done in the future. 

In summary, a semi-classical wobbling model is pre¬ 
sented about an extension of the original wobbling model 
by Bohr and Mottelson to derive the wobbling potential 
energy. This potential is anharmonic and described by 
a wobbling angle 0 , which is introduced as a dynamical 
variable to represent the semi-classical wobbling motion. 
The re-quantisation is made with this potential and the 
dynamical variable. The quantum equation for the an¬ 
harmonic wobbling motion is then derived. However, this 
anharmonic potential cannot reproduce the experimental 
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spectrum. In order to explain the discrepancy, a further 
anharmonic model taking into account a phase transition 
going into TAR from PAR is proposed as a possible in¬ 
fluences coming from the intrinsic microscopic degrees 
of freedom and the associated quantum tunnelling ef¬ 
fect. With the proper choice of the control parameter, 
it is demonstrated qualitatively that such a strong an- 
harmonicity seen in experiment can be well reproduced 
as a consequence of the spontaneous symmetry breaking. 
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